Trapped ion quantum computation with transverse phonon modes 
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We propose a scheme to implement quantum gates on any pair of trapped ions immersed in a 
large linear crystal, using interaction mediated by the transverse phonon modes. Compared with 
the conventional approaches based on the longitudinal phonon modes, this scheme is much less 
sensitive to ion heating and thermal motion outside of the Lamb-Dicke limit thanks to the stronger 
confinement in the transverse direction. The cost for such a gain is only a moderate increase of the 
laser power to achieve the same gate speed. We also show how to realize arbitrary-speed quantum 
gates with transverse phonon modes based on simple shaping of the laser pulses. 
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PACS numbers: 03.67.Lx, 32.80.Qk, 03.67.Pp 

Trapped ions have been demonstrated as one of the 
most promising systems for implementation of quantum 
computation. Different theoretical schemes have been 
proposed for quantum gate operations Q, 0> H El IS El > 
and many building blocks of quantum computing have 
been experimentally demonstrated H 0, 0, 0, 0, 
Il3l IhH fl5| . In an ion trap quantum computer, entan- 
gling gates between different ions are mediated through 
phonons in the collective ion motion. In all previous gate 
schemes H 0, H , the longitudinal phonon (LP) 
modes are exploited by kicking the ions along the axial 
direction of a linear trap. 

In this work, we propose to use the transverse phonon 
(TP) modes for gate operations. Compared with the 
conventional schemes (hereafter referred to as LP gates), 
gates involving TP modes (TP gates) have the following 
distinctive features: First, due to the strong confinement 
in the transverse direction, the TP gate is much less sensi- 
tive to ion heating and thermal motion. Even if the axial 
ion oscillation amplitude is significantly greater than the 
optical wavelength (outside of the Lamb-Dicke regime) , 
high-fidelity gates through the TP modes are still pos- 
sible. If (3 denotes the ratio of the center-of-mass (CM) 
trap frequencies for the transverse and the longitudinal 
directions (/? 3> 1 in typical experiments), we show that 
gate infidelity due to thermal ion motion is reduced by 
a factor ranging from (3 4 to /3 6 , depending on details of 
the initial ion temperature and the heating mechanism. 
This improvement may be particularly significant for a 
system of many ions or for ions confined in a microtrap 
[lfij . where ion heating and thermal motion may domi- 
nate gate erors. Second, the cost of using the TP modes 
is moderate, even though it is more difficult to excite the 
TP modes due to their strong confinement. For TP gates 
to have the same speed as LP gates, the intensity of the 
driving laser needs only to be increased by a factor of 
\J /3/2, a small factor when compared with the improve- 
ment in gate fidelity. Finally, we show that TP quantum 
gates can be operated with arbitrary speeds. Although 
the frequency splitting of the TP modes is significantly 
smaller than that of the LP modes, this does not impose 
any limit to the gate speed. High-fidelity fast TP gates 



are still possible through control of a simple sequence of 
laser pulses, which typically involves excitation of many 
TP modes. 

Transverse phonon (TP) modes. To design TP 
quantum gates, we first describe the structure of the TP 
modes. We consider a system of N ions confined in a 
linear trap. The phonon modes are obtained through 
diagonalization of the Hamiltonian for the ion exter- 
nal motion 17j. The ion motional Hamiltonian has 

the form H = £j=i |m ^ wnere M is the 
ion mass, and pt j is the momentum operator of jth 
ion along the direction £ (£ = x,y,z). The potential 
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counts for the Coulomb interaction between the ions as 
well as the external trapping potential, where r n j de- 
notes the distance between ions n and j, and is the 
CM trap frequency along the direction £. Typically, 
uj x ~ uj y 3> uj z , and one has a linear geometry with an 
ion chain along the z axis when tJ Xt y/u) z is larger than the 
critical ratio about 0.73-/V 086 18] . The equilibrium posi- 
tions z® (ar° = y° = 0) of the ions are determined by the 
condition [dV/dz n } Zn=z o = 0, (1 < n < N). With the 
dimensionless parameters u n = z^/£ (n = 1,2, ...,N), 
where I = y/ e 2 / '4tt€oMuj z sets the scale of ion spacings 
in the linear crystal, this condition yields a set of equa- 

„) 2 = 0, 



tions Uj-YfnJiViUj 



which can be solved numerically to determine z° for any 
larger ion array [l7| . To find all the phonon modes, 
one then just needs to expand the potential V around 
the ions' equilibrium positions ^ = z°) with 
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|. Under the harmonic approximation, 
we have V = (1/2) £ e , nJ dq^V/d^d^^o = 
n . A'l . where the matrix elements 
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(3, 



N 

E at 1 



~ a tl 



>? (n = j) 
(n + j), 
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with /?£ = lu^/lUz, a x — a v = —1, and a z = 2 
[Tij. The eigenfrequencies u^k = \/^tjk w t an( i e ig en_ 
vectors b|' fe of the normal phonon modes are obtained 
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from diagonalization of the matrix = [A^] with 
J2n A ij h n k = A «,fc b i' fc - Then, with the substitution 

= J2k h f k V^/ 2Muj ^A a i,k + 4,fc)> tnc motional 
Hamiltonian Hq reduces to the standard form Hq = 
SaLi ^ w ?,fe( a 5,fc a 4,fc + 1/ 2 )j expressed by the annihi- 
lation and creation operators a^,k> a \ k °f ^ ne ^h normal 
mode in the £ direction. 

In order to visualize the TP and LP modes, we plot 
the complete mode spectrum for a 10-ion array in Fig. 1 
(the modes along x and y directions are degenerate, so it 
is enough to show only the x- modes) . We choose the trap 
frequency ratio j3 x — 10, which is typical for experiments 
and larger than the critical value of 5.3 to stabilize a 
linear configuration for N — 10 ions. As opposed to the 
LP modes, the highest frequency TP mode is the center- 
of-mass mode at oj x . The frequency splitting between the 
CM mode and the second-to-highest mode (the bending 
mode) is about 0.05w z for N = 10, which is significantly 
smaller than the splitting (vo-l)u t of the corresponding 
LP modes (the spectral structure of the TP modes is 
inverted compared to the LP modes, as seen in Fig. 1). 
For entangling local ions (such as neighboring ions), it is 
best to use the low-frequency TP "zigzag" mode 20] as 
it is more resolved from the other TP modes and most 
insensitive to the ion heating. But the CM mode has 
the advantage that it is equally coupled to all the ions, 
and thus more appropriate for gates between nonlocal 
ions (such as ions at different edges of the chain) . When 
comparing features of gates using TP or LP modes, we 
parameterize the comparision of the CM modes for both 
cases. 
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FIG. 1: The spectrum of the longitudinal (LP) and transverse 
phonon (TP) modes for a 10-ion array. The dashed lines (No. 
1) denote the center of mass (CM) mode. 



In this case, the Hamiltonian for the laser-ion interaction 
has the form 

JV 

H = m 3 cos(Ak • q 3 + jj£)a? , (2) 

i=i 

where Ak • = Ak^g|(t), and flj denotes the two- 
photon Rabi frequency of the jth ion, which is propor- 
tional to the intensity of the driving laser. For conve- 
nience, fL is assumed to be real, but it can be time- 
dependent. 

Now we assume that the relative wave vector Ak is 
chosen along a certain direction £ (£ = x or z), and 
motion in all modes in this direction is in the Lamb- 
Dicke regime ?7f,fe\/%>fe + 1 *C 1 for all k, where = 
|Ak| ^Jhj2Muj^k is the Lamb-Dicke parameter and n^k 
the mean phonon occupation number of mode (£, k). 
Note that for TP quantum gates (£ = x), the lower fre- 
quency LP modes £ = z (as well as the other transverse 
mode £ = y) are decoupled and hence need not be con- 
fined within the Lamb-Dicke regime. To lowest order 
in rj^k and under the rotation-wave approximation, the 
interaction-picture Hamiltonian of the system is 

N 

j,k=i 

where the coupling constant g^j = ry^fcb|' , and Xj(t) = 
flj sin(/xt) is proportional to the state-dependent force on 
the j-th ion. 

The evolution operator corresponding to the Hamilto- 
nian Hi is given by 0, 0] 

U(t) = exp[* + 1 E 4( T )°K]> ( 4 ) 

3 3<n 

where the displacement operator 0j( T ) = 

T, k [ a lj( T ) a lk + »f*K T K. fc ] with a lj( T ) 

Jo Xj(t)9£ je luJ £- kt dt, and the conditional phase </>|„(t) = 

2/ r dt 2 J Q t2 dhJ^kXsi^gudinXnit^smuj^kih - h). 

In the case of a time-independent Qj , a* ■ (r) and </>j n (r) 
have the following explicit expressions: 



General formalism of trapped ion quantum o n k /,, , „iu e nrt ,, „—,/-, ,_\ 1 „•, , cW.^U 

gates, first, we give a general tormalism tor multi-ion a f = ■ 

entangling quantum gates, which is valid with either TP 
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or LP modes. The qubit for each ion is represented by 
two hyperfine states, denoted as |0) and |1) in general. gk , , 2 it + sin(2 r)l 

The gate is achieved b y apply ing a state-dependent ac- t/>$ n — 20jfi„ 2 ^' 3 ' ^ < ^ — — — — ^— 

Stark shift on the ions [lll[l5|, using two laser beams of fe " £>k L " 

equal intensity, wavevector difference Ak and frequency 
difference fi. As it is common in experiments, we assume 

that the average ac-Stark shifts are the same for the |0) C0S H,fcT) smprr) - ^cos(^t) sin(^ ;fc r)] I ^ 

and |1) states for the ions in their equilibrium positions. ~ ^£,k J 



(6) 
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To drive a conditional phase flip (CPF) gate between 
arbitrary ions j and n, we take Qj to be nonzero only for 
these two ions, and set the laser detuning [i and the gate 
time r so that </>|(r) = <^4(r) = and (j£j n (j) = ir/4. 
In this case, the evolution operator U(t) reduces to the 
CPF operator Uj n = exp(«7rcr|(T^/4). 

TP versus LP quantum gates. To compare quan- 
tum gates based on TP vs. LP modes, let us start with 
the assumption that the driving laser can address individ- 
ual phonon modes through frequency selection (resolved- 
sideband addressing). This requires the gate time r to 
be much larger than ta = 27r/A, where A is the charac- 
teristic frequency splitting of the phonon modes (we will 
see later that only one phonon mode dominates when 
t > 2ta). The sideband addressing assumption, al- 
though not essential, allows us to derive simple analytic 
relations that permit a direct comparison between TP 
and LP gates. 

With sideband addressing, we dominantly excite a par- 
ticular phonon mode (£,p) with frequency uj^ tP by ad- 
justing the laser detuning /i close to with fi = 
lj^^p + 2ttI^/t, where is an integer, typically chosen 
as 1 or — 1 0, ■ The qubit state and the motion state 
should be disentangled at the end of the gate, which re- 
quires </>|(r) = </4(t) = and thus r = V^k/lj^^ (see 
Eq. (4)), where li is another integer (typically, 11 ^> I 



From Eq.@), the conditional phase shift is found to be 

(7) 
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The condition 4>j n ( T ) = can be satisfied with an 
appropriate choice of the Rabi frequency O. 

To consider inherent infidelity of the gate operation, 
we note that all the above results are derived based on 
the Lamb-Dicke condition. In practice, the Lamb-Dicke 
parameter is finite, and the thermal motion of the ions 
will induce some fluctuation of the Rabi frequency il 2 in 
Eq. (6) and lead to gate errors. To estimate this noise, 
we need to expand the Hamiltonian of Eq. J5J to higher- 
orders in the Lamb-Dicke parameters. The effect of these 
higher-order terms is to replace f2 in Eq. (6) with an 
effective Rabi frequency that depends on the phonon 
number of the ion vibrational mode (£,p). To the next 
order of the Lamb-Dicke parameter 77{ lP , we find (fi^) 2 ~ 
Of [1 - r?f, p (2n € + 1)] 3]. The gate fidelity Fj |2 is then 

found to be F% = 



Kp(2"+1) 



g 2^2 Z^n=0 

When the initial phonon number distribution Pj takes 
the form of a thermal state, Pj = n^/(l + n^) n+1 , we 

find to lowest order in r)£ tP a gate infidelity Ff n = 1 — P| 
of 



1/8). 



(8) 



t ? -r 

As the TP mode has a larger vibrational frequency, the 
TP quantum gates have a significantly smaller gate infi- 
delity from thermal ion motion. Even if we assume the 



TP and the LP modes have the same mean phonon num- 
ber n x = n z , the infidelity for the TP gate is smaller by 
a factor of Ff n /F? n = T)£, p /vi,p = (wz, P /wx,p) 2 - In prac- 
tice, if the TP and LP modes are subject to the same 
heating mechanism, and initially prepared with the same 
laser cooling technique, we further expect n x <C n z due 
to differences in the initial temperature and the heat- 
ing rate of TP vs. LP modes. The temperature limit 
Tl can be considered as independent of the phonon fre- 
quency uj^^p for Doppler cooling, and is roughly propor- 
tional to l/w{, P for the Raman sideband cooling. So the 
contribution of Tl to the mean phonon number nj, es- 
timated as kBT^/hui^^, is taken to be (1/u^p) 7 , where 
7 is between 1 and 2. The ion heating rate for the 
phonon mode (£,p) is proportional to the noise power 
spectrum S(u>^ tP ) at the frequency w^ iP |24| . taken to be 
independent of frequency (white noise) or proportional 
to l/w£, P (1// noise). For these practical noise sources 
|24j |. the average phonon occupation number therefore 
scales as (l/^ iP ) 7 , with 7' again between 1 and 2. If 
we assume the term n 2 dominates in the infidelity ex- 
pression (8), which is likely for many ions in a crys- 
tal, the infidelity ratio of TP vs. LP gates is then 

F fn/ F in ~ nlvij,/ (nlvip) ~ (^, P M,p) 2+7+7 ', where 
we have assumed that the gate time r is the same for 
both cases. For the CM modes, co z ,p/^x,p is given by the 
trap frequency ratio f3 x = ui z /w x . So, compared with LP 
gate, the inherent infidelity of the TP gate could be re- 
duced by a factor of (3 X to which is significant even 
for a moderate trap frequency ratio of (3 X ~ 5. 

Now we look at the cost of the TP quantum gate. As 
the TP modes have a higher vibrational frequency, it is 
harder to excite them, and we need more laser intensity 
for the same gate speed. From Eq. (7), to have the same 
gate time r, the ratio of the required laser intensity I x /I z 
(note that 1^ is proportional to the two-photon Rabi fre- 
quency fi^) is given by I x /I z ~ T] z , p /r) XtP = ^u x ,p/u) z , P = 
y/JJ^ (we neglect l^/l'^ in Eq. (7) as it is typically much 
less than 1). So, although we need additional laser power 
for the TP quantum gate, this cost is moderate compared 
with the improvement we achieve in the gate fidelity. If 
we take into account of different laser excitation config- 
urations for the TP and the LP gates, this cost is even 
less. For the LP gate, the relative wave vector Ak of 
the Raman laser beams needs to be along the trap axis, 
but one cannot directly apply lasers in that direction, so 
in practice both of the laser beams have a 45° to the 
trap axis 0,0]. However, for the TP quantum gate, 
one can apply counter-propagating laser beams along the 
x-axis so that Ak is perpendicular to the ion string. Be- 
cause of this difference, the laser intensity ratio I x /I z is 
actually ^ f} x /2 instead of V/^r- The above laser con- 
figuration difference also gives some practical advantages 
for the TP quantum gates: first, as the laser beams are 
perpendicular to the ion string, it is easy to have the 
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same relative laser phase for different ions, as assumed 
in the gates above. For LP gates, it is difficult to achieve 
such a condition for different ions as they are not equally 
spaced. It usually requires subtle control of the ion dis- 
tance through adjustment of the trap frequency [ToL fTl| ■ 
and such a technique is not scalable to many ions. Sec- 
ond, with the transverse focused laser beams, it is also 
easier to achieve separate addressing of different ions. For 
the TP gate, one does not need to have a large longitu- 
dinal uj z to achieve the Lamb-Dicke condition, so the ion 
distance is not limited, and one can combine separate 
addressing with a many-ion setup, which is desired for 
scalable quantum computation. 

Arbitrary-speed TP gates through minimal 
control of the laser beams. As we have mentioned 
before, the TP modes have small frequency splittings, 
so resolving a particular TP mode could be difficult for 
a large ion array. If we want to achieve a high-speed 
gate, it is necessary to go beyond the sideband address- 
ing (single- mode) limit. Fortunately, for any practical 
qubit number (up to a few hundreds, for instance), it is 
always possible to take into account of all the phonon 
modes and design high-fidelity gates with no limitation 
to the gate speed 0, IE • These fast gates are in general 
based on control of the laser pulses. 

Here, similar to Ref. Q, we use a simple sequence 
of laser pulses which take minimum experimental con- 
trol. Wc chop a continuous-wave laser beam into m 
equal-time segments (through acoustic/electric modula- 
tor, for instance), with a constant but controllable Rabi 
frequency £l p for the pth (p = I, 2, • • • , to) segment. The 
state-dependent force in the Hamiltonian (3) then 
takes the form x(t) = £! p sm(fit) for the time interval 
(jp — \)t / m < t < pr/m (For simplicity of the notation, 
we omit the direction index £ in the following and take 
£ = x by default). With a large number of ions but 
a small number of control parameters f2 p , the displace- 
ments otj{r) (and thus 4>j{r)) in the evolution operator 
(4) may not exactly reduce to zero. But as long as they 
are small, we still can get a high-fidelity gate. In this 
case, for a gate on the ion-pair (J, n) in an TV-ion array, 
the gate infidelity due to the residue small displacements 
a* is given by 

i^«£^[K fc | 2 + k;| 2 ]/4, (9) 

k 

where j3k = cothj-y/ fik /4Zn(l + l/fix)] with fi\ repre- 
senting the mean phonon number of the CM TP mode. 
For the above expression, we have assumed that the TP 
modes are in their thermal equilibrium state. The task 
of control is to get a small infidelity F„ by choosing an 
optimal laser detuning \i and a minimum number of the 
control parameters Q p . 

First, let us consider the gate with a single continuous- 
wave laser beam (the number of segments to = 1). The 




FIG. 2: For the two edge ions in a 10-ion array, the TP gate 
fidelity as a function of the detuning fi with the gate time 
r = 5to, f%i = 3, f3 x = 10, and the number of segments 
m — 1,3, 5, 8, respectively. 



only control parameter is the detuning /i. We find that 
for the gate on two edge ions (1st and 2nd ions) in a 10- 
ion array, as long as the gate time r > 37tq, where tq = 
2ir/ui z , the gate infidelity Fi n < 0.99%. For this and the 
following calculations, we take fi\ ~ 3 which corresponds 
to a pretty high temperature. The optimal /j, is very close 
to the value ui x + 2it/t. Note that a gate with r « 37ro 
has been faster than any ion gate implemented so far in 
the lab |ll|,|l5|. For this gate, the time r is close to ta = 
27r/A (r ps 1-8ta), where A is the frequency splitting 
between the CM TP mode and the bending mode. So, 
besides the dominant CM mode, various TP modes are 
indeed slightly excited during the gate and contribute to 
the conditional phase ^>j n (r). But with the optimal pL, 
all these modes evolve along an almost-closed loop in the 
phase space (cq w 0, although not exactly), so we still 
have a high fidelity gate. The required Rabi frequency 
from the exact numerical calculation is very close to the 
one given by the analytic formula (7) from the single- 
mode approximation (within a 6% error). 

If we further increase the gate speed with r < 37to, 
the gate fidelity quickly decreases, so we need to chop 
the continuous-wave laser beam into more segments with 
to > 1 to increase the fidelity. With a sufficiently large 
to, a high-fidelity gate with an arbitrary gate speed 
can be achieved. In Fig. 2, we show the calculation re- 
sult for the gate time r = 5tq. With the number of 
segments to = 1,3,5,8, the gate infidelity is given by 
10%, 4.9%, 1.0%, 0.1%, respectively, with the optimized 
parameters /i and - ■ ■ , Q m . We also calculate 

the gate infidelity for other ion-pairs, and the results 
are qualitatively similar. For instance, with t = 5to 
and the number of segments to = 1,3,5,8, respectively, 
the gate infidelity F in is given by 5.5%, 1.8%, 0.22%, 
0.07% for the two center ions (5th and 6th ions), and by 
40%, 25%, 8.5%, 0.99% for the 1st and 10th ions at the 
far ends of the string (the worst case). 

In summary, we have proposed to use the transverse 
phonon mode to achieve quantum gates, and shown that 
these TP gates are much less sensitive to ion heating and 
thermal motion compared with conventional gates based 
on the longitudinal modes. We also show how to achieve 
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arbitrary-speed TP quantum gates based on control of 
the laser pulses. 
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